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ABSTRACT 

We use the cannonball model of gamma ray bursts (GRBs) and public data 
from the first day of observations of GRB 090618 to predict its X-ray and optical 
lightcurves until very late times, and, in particular, the emergence of a photo- 
metric and spectroscopic signature of an SN akin to SN1998bw in its optical 
afterglow with an anticipated peak brightness of magnitude ~ 23.2 in the R band 
around July 10, 2009, if extinction in the host galaxy can be neglected. 



1. Introduction 

The relatively nearby bright gamma ray burst (GRB) 090618 that was discovered by the 
Swift Broad Alert Telescope (BAT) on June 16, 2009 at 08:28:29 UT (Schady et al. 2009a) 
provides another good opportunity to investigate the association of long duration GRBs 
with supernova explosions and to test theoretical models of long GRBs (e.g. Malesani 2009, 
Dado and Dar 2009). Two such models have been used extensively to analyze GRBs and 
their afterglows (AGs), the fireball (FB) model (for recent reviews see, e.g., Meszaros 2006; 
Zhang 2007) and the cannonball (CB) model (e.g. Dado & Dar & De Riijula 2002, hereafter 
DDD002, Dar & De Rujula 2004, hereafter DD2004; Dado & Dar & De Riijula 2009a,b, 
hereafter DDD2009a,b, and references therein). The two models are quite different and at 
most only one can correctly describe GRBs. Until recently, the FB model has been widely 
accepted as that one. However, the rich data on GRBs and their afterglow accumulated 
from space based observations, in particular from the Swift and Fermi satellites, from early 
time observations with ground based robotic telescopes and from late-time follow-up obser- 
vations with large telescopes, have challenged this prevailing view (DDD2009a and refer- 
ences therein): Synchrotron radiation (SR), the only radiation mechanism in the original FB 
model, cannot explain simultaneously the prompt optical emission and the hard X-ray and 
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gamma-ray emission from GRBs which were both well measured in some bright GRBs such 
as 990123 and 080319B. The prompt hard X-ray and gamma-ray pulses cannot be explained 
by synchrotron radiation from internal shocks generated by collisions between conical shells. 
Neither can SR explain their typical energy, spectrum, spectral evolution, pulse-shape, rapid 
spectral softening during their fast decay phase and the established correlations between 
various observables. Moreover, the early time high energy emission is uncorrelated to the 
prompt keV-MeV emission and lags behind it. As for GRB afterglow (AG), contrary to the 
predictions of the FB model, the broad band AG of GRBs is highly chromatic at early times, 
the AGs of the brightest GRBs do not show jet breaks, and, in canonical AGs where breaks 
are present, they are usually chromatic and do not satisfy the closure relations expected from 
FB model jet breaks. 

The situation concerning the CB model is different. The predictions of the model 
which were derived in fair approximations from its underlying assumptions were shown to 
describe correctly the main observed properties of GRBs and to reproduce successfully the 
diverse broadband lightcurves of large representative sets of both long GRBs and XRFs 
(e.g., DDD2009a and references therein) and short hard bursts (DDD2009b), in particular 
of relatively nearby bright GRBs with precise and well sampled lightcurves. Here we repeat 
the exercise of predicting from early time observations the entire late time behaviour of the 
afterglow of the recent, relatively nearby, bright long GRB 090618 and the emergence of a 
photometric and spectroscopic signature of an SN akin to SN1998bw in its lightcurve before 
their observations. 



2. GRB090618 

At 08:28:29 UT, the Swift Burst Alert Telescope (BAT) triggered and located the bright 
long duration gamma ray burst (GRB) 090618 (Schady et al. 2009a) at redshift z=0.5A 
(Cenko et al. 2009b). About 90% of the GRB energy measured by BAT was emitted within 
T90=113s (Baumgartner et al. 2009). The Swift X-ray telescope (XRT) began follow up 
observations of its X-ray lightcurve (see Fig. 1) 124 s after the BAT trigger and its UVO 
telescope detected its optical AG 129s after trigger (Schady et al. 2009b ). The burst was 
detected also by AGILE (Longo et al. 2009), Fermi GBM (McBreen et al. 2009), Suzaku 
WAM (Kono et al. 2009), KONUS-WIND and KONUS-RF (Golenetskii et al. 2009). The 
burst light curve showed a smooth multipeak structure with 4 prominent peaks (one followed 
by three much brighter overlapping peaks) with a total duration of 160 s. Significant 
spectral evolution was observed during the burst. The spectrum at the maximum count rate, 
measured from T+62.720 to T+64.0 s, was well fitted (Golenetskn et al. 2009) in the 20 keV 
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- 2 MeV range by the Band function (Band et al. 1993) with a low-energy photon index -0.99 
(-0.06, +0.07), a high energy photon index -2.29(-0.5, +0.23), and peak energy Ep=440t70 
keV while the time integrated spectrum had a low-energy photon index -1.28 ±.02 , a high 
energy photon index -2.66(-0.2, +0.14), and a peak energy Ep=186il:8 kcV. The isotropic 
equivalent energy in the 8-1000 keV band was Eigg — 2.0 x 10^^ erg (standard cosmology). 

The bright optical afterglow of GRB090618 was first detected by the ROTSE III robotic 
telescope 23.9 s after the BAT trigger (Rujopakarn et al. 2009) and by the Palomar 60- 
inch telescope (Cenko ct al. 2009a), the Katzman Automatic Imaging Telescope (Perley ct 
al. 2009, Li et al. 2009) and the UVO Telescope aboard Swift (Schady et al. 2009a) within 2 
minutes after trigger. Absorption features which were detected in its bright optical AG with 
the 3m Shajn telescope at Lick observatory yielded a redshift of z=0.5A. Its optical AG was 
followed up by many telescope and reported shortly after in GCN circulars (sec Table I). Its 
R band light curve reported in these GCN circulars before June 27, 2009 is shown in Fig. 2 . 



3. The CB model 

In the cannonball (CB) model (DDD2002, DD2004, DDD2009a, and references therein) 
long-duration GRBs and their AGs are produced by bipolar jets of highly relativistic CBs of 
ordinary matter which are ejected (Shaviv & Dar 1995, Dar & Plaga 1999) in core-collapse 
supernova (SN) explosions akin to SN 1998bw (Galama et al. 1998). Their prompt MeV 
gamma-rays and hard X-rays are produced by the thermal electrons in the CBs' plasma 
by inverse Compton scattering (ICS) of glory photons - photons emitted/scattered into a 
cavity created by the wind/cjccta blown from the progenitor star prior to the SN. When 
the CBs cross the wind/ejecta and coast through the interstellar medium (ISM) behind it, 
the electrons of the ionized gas in front of them that arc swept in and Fermi accelerated 
by the CBs' turbulent magnetic fields emit synchrotron radiation (SR) which dominates the 
'prompt' optical emission and the broad band afterglow emission. ICS of the SR radiation 
by these electrons and the decay of 7r°'s produced in collision between the swept-in wind 
and ISM protons and the ambient CB protons produce the 'prompt' high energy emission 
simultaneously with the optical emission. Within the CB model, the above radiation mech- 
anisms with the burst environment suffice to provide a sufficiently accurate description of 
the observed radiations from GRBs at all times and all detected wavelengths. 
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3.1. The optical light curve 



In the CB model, the observed optical light has three origins: the ejected CBs, the 
SN explosion, and the host galaxy. The optical light curve is the sum of their energy flux 
density: 

Fag [i^, t] = Fob [v, t] + Fsn [y, t]+FHG[i^A- ( 1 ) 

The contribution of the host galaxy, F^g^ is usually extracted from very late time observa- 
tions when the CB and SN contributions become negligible. In the case of GRB090618, a 
faint object is visible in the SDSS r and i frames at the position of the optical afterglow, 
likely its host galaxy. Compared to nearby SDSS stars, its r magnitude was estimated by 
Malesani (2009) to be 22.7±0.3. 

The energy flux density of an SN like SN1998bw with an energy flux density -Fbu,[z/, t] at 
redshift 2;;,^ = 0.0085 (Galama et al. 1998) placed at a redshift z is given by (e.g., DDD2002), 



FsNii'.t] = k 



Dliz) 



Fbw[kv-, t/k] , 



(2) 



where k — {l-\-z)/{l-\-Zbw), -^{1^,^) is the attenuation of the observed SN light at frequency 
V arriving along the line of sight, and Dl{z) is the luminosity distance to redshift z (we use 
the standard cosmology with Qm = 0.27, Qa — 0.73 and ifo = 71 kms~^ Mpc~^). 

In the CB model (DDD2009a,b and references therein) the afterglow emission begins 
when the CBs encounter the wind/ejecta of the progenitor star. It is dominated by syn- 
chrotron radiation (SR) from the electrons of the ionized wind and interstellar-medium in 
front of the CBs which arc swept in and Fermi accelerated by the CBs' turbulent magnetic 
fields which we assume to be in approximate energy equipartition with their energy. 

The SR, isotropic in the CB's rest frame, has a characteristic frequency, Ub{t): the typical 
frequency radiated by the electrons that enter a CB at time t with a relative Lorentz factor 
7(i), and spiral around its equipartition magnetic field (with the incident protons). In the 
observer's frame: 
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1/2. 



Hz, 



(3) 



1 + Z 1012 

where uq ~ 3.85 x 10^^ Hz ~ 160 eV and S{t) is the Doppler factor of the CB. The spec- 
tral energy density of the SR from a single CB at a luminosity distance Dl is given by 
(DDD2009a): 
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where p~2.2 is the typical spectral index of the Fermi accelerated electrons, and A(z/, t) is 
the attenuation of photons of observed frequency v along the line of sight through the CB, 
the host galaxy (HG), the intergalactic medium (IGM) and the Milky Way (MW): 

A{v,t) = exp[-r,(CB)-r,(HG)-r,(IGM)-r,(MW)]. (5) 

The opacity t^{CB) at very early times, during the fast-expansion phase of the CB, may 
strongly depend on time and frequency. The opacity of the circumburst medium [r,^(HG) 
at early times] is affected by the GRB and could also be t- and i/-dependent. The opacities 
ri/(HG) and T,y{lGM) should be functions of t and u, for the line of sight to the CBs varies 
during the AG observations, due to the hyperluminal motion of CBs. 

3.2. The early-time SR 

At early-time, before the CB has swept a mass comparable to its rest mass both 7 and 
6 stay put at their initial values 70 and 6o- Then, Eq. (jll) yields an early-time SR light 
curve, FsR[i',t\ (xe~'^w R'^n^^^^^^'^ v"^ . Since rcxt, a CB ejected into a windy density profile, 
n oc created by the mass ejection from the progenitor star prior to its SN explosion, 
emits SR with an early-time light curve, 

Fsr[^, t] oc '-^^ z.-^ ^r(^^) u-^ . (6) 

For a CB ejected at time t,, the time t must be replaced by t — tj, the time after ejection. 
In the 7-ray and X-ray bands, the SR emission from a CB is usually hidden under the 
prompt keV-MeV ICS emission. But, in many GRBs, the asymptotic exponential decline 
of the energy flux density of the prompt keV-MeV ICS emission is taken over by the slower 
power-law decline, FsR[i'.,t](xt~^^ with rx=/3x+l~2.1 of the synchrotron emission 

in the windy ~l/r^ circumburst density before the CB reach the constant ISM density and 
the AG enters a plateau phase (see examples, e.g., in DDD2009a). Note that the 'prompt' 
optical emission that is dominated by SR, decays initially like FsR[i',t] cxt^^'^ u^-^ since the 
spectral index of the unabsorbed SR emission at frequencies well below the bend frequency 
is /5o~0.5 . As the wind density decreases with distance, the bend frequency may cross the 
optical band while the CB is still in the wind, yielding a steeper decay, FsrIv, t] -^f"^-^ u^^-^ . 

3.3. The plateau, the break and the late time decay 

During the coasting phase of a CB in a constant density ISM the behaviour of its SR 
lightcurve as given by Eq. (jlj) is dominated by the time dependence of its Lorentz factor 
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7(it). From energy-momentum conservation it follows that 



7(t) 



7o 



with 



(7) 
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This deceleration law is for the case in which the ISM particles re-emitted fast by the CB 
are a small fraction of the flux of the intercepted ones. As can be seen from Eq. ([Tj), 7 and 
5 change little as long as t <^ t;, = [1+ 7o 6*^]^ to which results in the shallow decline/plateau 
phase of the AG. In terms of typical CB-model values of 70, -R, and n, 
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For t^tb, 7 and S decrease like t"^/^ . The transition 'j(t) ~ 70 — 7 ~ 70 {t/to)~^^^ induces 
a bend (the so called 'jet break') in the synchrotron AG from a plateau to an asymptotic 
power-law decay, 



Fsr[u, t] OC t-P/2-V2 ^ ^-/3 ^ ^-r+1/2 ^-T+l ^ 



(10) 



with a power-law in time steeper by half a unit than that in frequency. 



4. The X-ray lightcurve of the afterglow of GRB090618 

The X-ray lightcurve of GRB090618 (Evans et al. 2008) shows the canonical behaviour 
predicted by the CB model (e.g. DDD2002; DDD2009a) and displayed by many Swift GRBs 
(e.g., Nousek et al. 2006). This behaviour is well reproduced by the CB model as shown in 
Fig. 1. The fast decline phase with a rapid spectral softening is that predicted for ICS of 
glory light (e.g. Dado, Dar and De Riijula 2008a; DDD2009a): FicsW,t]^t-'^ e'^^^/^^^p. 
where t and tp are measured relative to the beginning of the last large prompt emission 
episode (we used a best fit value, tp = 12.9). The sharp transition around 300 s to a shallow 
decline/plateau phase with a constant hardness ratio is produced when the synchrotron 
afterglow given by Eq. (jlj) takes over. The shape of the SR afterglow of GRB090618 was 
reproduced with three best fitted parameters, 70 ^ = 1.10 and to = 312s which yield = 1540s 
and p = 2.08, This best fit value of p satisfies well the CB model relation, p = 2r— 2 = 2.02ild0.10, 
where we used the photon spectral index r=2.008 (+0.047, -0.046) that was reported in the 
Swift repository (Evans et al. 2008). The 'jet break' takes place when the jet of CBs gathers 
a mass comparable to its rest mass (DDD2002; DDD2008b; DDD2009a) as given by Eqs. ([8]) 
and Q). The post-break power- law decay of the AG is well described (DDD2008b) by 
Eq. (do]). 
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5. The optical AG of GRB 090618 and emergence of an underlying SN? 

The late time optical AG can be estimated by extrapolating the post break power- 
law decay to the time of the anticipated emergence of an SN signature akin to that of 
SN1998bw at the burst location. However, unlike the X-ray light curve which was inferred 
from a continuous follow-up measurements with the same telescope (Swift XRT), the optical 
lightcurve of GRB090618 constructed from reported measurements in GCNs with different 
telescopes at different times, locations, atmospheric and seeing conditions, calibrations and 
spatial resolutions. In particular, the detection of the SN signature depends on a precise 
subtraction of the host galaxy contribution to the observed lightcurve. In view of all that 
we preferred to best fit the early-time observational data on the optical afterglow which 
were reported in the GCNs listed in Table I and to use the CB model with the parameters 
determined from the X-ray AG and the early time observational data (Table I) on the optical 
AG of GRB090618 to predict the late time R-band lightcurve of the optical transient. A 
host galaxy contribution of r=22.7±0.3 (Malesani et al. 2009) was subtracted from the last 
two data points, and the anticipated SN1998bw-like contribution at the host location was 
dimmed by the Galactic extinction along the line of sight corresponding to E(B-V)=0.09. 
The results of this exercise are presented in Fig. 2. 

Although the very early optical emission does not directly affect the late-time behaviour 
of the AG, we have also fitted the very early ('prompt') optical emission in a windy environ- 
ment as given by Eq. (Q, for completeness and in order to demonstrate the validity of the 
CB model. The fit to the two prompt emission peaks was obtained with /5o=0.50, tj=25.1s, 
a=1.47s, tea;p=19.4s for the first peak and tj=52.4s, a=98s, tea:p=40.3s for the second. The 
scarcity of data points during the rise part of the peaks allows for other equally good fits. 

A careful subtraction of the host galaxy contribution from the optical lightcurve of 
GRB090618 measured with large telescopes of good spatial resolution should have shown 
already (on June 29, 2009) spectroscopic and photometric evidence for an underlying super- 
nova if its brightness is similar to that of SN1998bw. As of today, the R band lightcurve 
should show a plateau with a very small rise towards a shallow maximum around July 10, 
when an SN akin to SN1998bw with only little extinction in the host galaxy will reach its 
peak brightness of magnitude ~ 23.2 in the R band. The detailed data on the broadband 
afterglow of GRB090618, which is less bright than that of GRB030329, can still provide 
another useful test of the long-duration GRB-SN association and of GRB models. 
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Table 1. GCNs Reporting Follow-up R-band Optical Observations of GRB090618 Before 

June 25, 2009 



GCN 


Telescope 


Reported by 


Comments 


9512 


Swift UVOT 


Schady et al. 


Finding exposure 


9513 


Palomar 60-incli 


Cenko et al. 




9514 


KAIT 


Per ley et al. 




9515 


ROTSE-III 


Rujopakarn et al. 


Earliest Detection 


9517 


KAIT 


Li et al. 


Follow-up 


9518 


Lick 3-m 


Cenko et al. 


Redshift Measurement 


9519 


OAGH 2.1-m 


Carraminana et al. 


Follow-up 


9520 


Faulkes 


Melandri et al. 


Follow-up 


9522 


Mt. Lemmon 1-m 


Im. ct al. 


Follow-up 


9526 


SDSS 


Malcsani 


Host Detection 


9531 


Liverpool 


Cano et al. 


Follow-up 


9539 


Shajn 


Rumyantsev et al. 


Follow-up 


9542 


SAO RAS 6-m 


FatkhuUin et al. 


Spectroscopy 


9548 


RTT 1.5-m 


Galeev et al. 


Follow-up 


9575 


SARA 0.9-m 


Updike et al. 


Follow-up 


9576 


Himalayan 2-m 


Anupama et al. 


Follow-up 
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Fig. 1.— Comparison between the Swift XRT lightcurve of GRB090618 (Evans et al. 2009) 
and its CB model description (see text). The time of the deceleration bend/break tb is 
indicated by a down-pointing arrow. 
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Fig. 2. — Comparison between the R band lightcurve of GRB090618 extracted from the 
GCN circulars hsted in Table I and its CB model description in terms of the parameters 
which were extracted from the CB model fit to its early X-ray lightcurve measured by the 
Swift XRT (Evans et al. 2009). The plotted SN hghtcurve is that of SN1998bw at the burst 
location. Its maximum brightness is expected around July 10. 



